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ABSTRACT 


In  this  research,  a  sequence  of  models  is  constructed  to  simulate  the  movement  of  fluids  and 
chemicals  in  the  cerebrovascular  system.  One  model  simulates  the  nonsteady  response  of  perfu¬ 
sion  in  various  sections  of  the  brain.  In  a  second  model,  certain  relevant  parts  of  the  body  are 
added  to  form  a  single  brain-body  model.  A  third  model  simulates  the  transport  of  selected  chemi¬ 
cal  components  through  the  cerebrovascular  system.  Predictions  derived  from  the  brain  model 
were  shown  to  be  well  within  the  range  of  available  clinical  observations.  The  brain-body  model 
describes  the  interaction  between  the  cerebral,  the  cardiovascular  and  the  respiration  systems.  It 
is  excited  by  expiration/inspiration  fluxes  and  accounts  for  the  effects  of  hydrostatic,  environmental 
pressures,  flight  maneuvers  with  excessive  (head  to  bottom)  gravitation  acceleration  and  resusci¬ 
tation  procedures.  In  simulating  chemical  processes  in  the  brain,  the  model  accounts  for  C02, 
HCOg'  and  H+,  as  they  are  transported  by  perfusion  and  diffusion  in  the  presence  of  chemical 
reactions.  This  model  also  focuses  on  the  flow  control  between  brain  arteries  and  capillaries,  due 
to  changes  in  C02  concentration. 


NON-STEADY  COMPARTMENTAL  MODEL  OF  INTERACTIVE  PERFUSION 
BETWEEN  CEREBRAL  AND  BODY  SYSTEMS 

S.  Sorek,1  K.  Allen,2  M.  Feinsod3  S.  Ben  Haim4 
J.  Beat5  and  L.  Bunt6 


ABSTRACT 

A  lumped  parameter  model  is  developed  to  simulate  perlusive  flux  and  pressure  interaction 
between  the  cerebral  and  the  body  systems.  Its  objective  is  to  study  the  dynamic  interaction 
between  the  cranio-spinal,  body  respiratory  and  the  heart  systems  that  influence  the  brain,  as  well 
as  the  influence  of  conditions  in  the  environment  on  the  body.  By  providing  forecasts  ot  depar¬ 
tures  from  a  normal  behavior,  the  model  will  serve  the  following  medical  purposes,  (a)  facilitate 
the  understanding  of  the  physiology  and  the  mechanisms  that  preserve  the  delicate  brain  in  the 
face  of  living  stress,  and  (b)  provide  information  for  management  In  deviant  cases. 

The  compart  mental  model  consists  ot  six  compartments  that  describe  the  cerebral  system, 
and  eight  compartments  that  are  assigned  to  the  body  system.  The  model  also  accounts  for  the 
surrounding  environment  affecting  the  abdominal,  inhale/exhale  fluxes.  Altitudes  assigned  to  the 
compartments  introduce  the  effect  of  hydroslatic  pressure. 

Key  Words:  Compartmental  modelling,  interaction  between  cerebral  respiratory  and  heart  sys¬ 
tems;  perfusion  ftux  and  pressure;  conductances;  compliances;  environmental 
pressure;  hydrostatic  pressure;  inhale/exhale  fluxes. 
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1.  INTRODUCTION 

In  earlier  papers  (Kami  at  al.  1987;  Sorek  et  al.  1987a,  1987b),  a  compart  mental  model  was 
built  to  simulate  the  perfusion  In  the  cerebrovascular  system.  The  model  provided  information  on 
pressures  and  fluxes  in  response  to  excitation  in  the  form  of  temporal  flux  and  pressure  changes 
at  the  internal  carotid  artery.  The  influx  to  and  efflux  from  the  cerebrovascular  system  were  intro¬ 
duced,  respectively,  as  external  conditions  imposed  on  the  internal  carotid  artery,  and  the  jugular 
bulb. 

In  the  present  work,  this  model  is  extended  to  include  both  the  cerebrovascular  and  the  body 
system  In  a  single  model.  The  objective  is  to  enable  the  study  of  the  dynamic  interaction  between 
these  two  systems.  In  this  model,  relevant  parts  of  the  brain  and  of  the  human  body  are 
represented  as  compartments  that  interact  with  each  other,  e.g.,  in  the  form  of  pressure  transmis¬ 
sion  and  exchange  of  fluid.  Each  such  compartment  is  represented  by  lumped,  or  averaged,  pro¬ 
perties  and  state  variables  of  that  organ.  Examples  of  such  lumped  parameters  are  pressure  in  a 
compartment,  influx/efflux  through  its  boundaries  and  fluid  source/sink  terms. 

The  compartmental  model  is  comprised  of  six  compartments  that  represent  the  cerebral  sys¬ 
tem  and  eight  compartments  assigned  to  the  body  system.  The  body  portion  of  the  compart  men¬ 
tal  design  involves  only  those  parts  of  the  body  that  are  relevant  to  the  cerebral  system,  i.e., 
interact  grossly  with  the  latter.  Accordingly,  these  parts  include  mainly  the  heart  and  the  respira¬ 
tory  systems. 

The  model  also  accounts  for  the  influence  of  the  environment  surrounding  the  body,  e.g., 
atmospheric  pressure,  pressure  at  high  and  low  altitudes,  high  underwater  pressure,  vacuum  and 
excessive  gravity  acceleration,  such  as  encountered  in  flight  maneuvers.  It  also  takes  into  con¬ 
sideration  differences  in  elevation  between  compartments.  This  means  that  H  takes  into  account 
hydrostatic  pressure. 

Excitations  of  the  model  are  introduced  through  changes  in  inhale/exhale  rates  and  environ¬ 
mental  pressure  changes  exerted  on  the  abdominal.  Extreme  intervention,  such  as  resuscitation 
and/or  clogging  In  the  heart  system,  can  be  introduced  via  changes  in  the  appropriate  compart¬ 
mental  conductances  and  compliances. 
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The  diagrammatic  layout  o(  the  compartments  is  such  as  to  assist  a  medical  clinician  in  the 
interpretation  of  the  model  image  of  the  human  body,  as  when  examining  an  X-ray  plate.  The 
model  wM  guide  clinicians  by  predicting  the  brain-body  perfusion  responses  to  various  excitation, 
thus  enhancing  the  physiological  understanding,  faut  finding  and  consequent  management  in 
deviant  cases  of  this  complex  system. 

2.  MODEL  TOPOLOGY 

The  cerebral  section  of  the  model  is  based  on  the  works  of  Kami  et  at.  (1987)  and  Sorek  et 
al.  (1987),  and  is  comprised  of  the  following  compartments  (Fig.  1): 

Arterial  Cranium  (Aq  )  -  Consists  of  four  supply  vessels  through  the  right  and  left  internal 
carotid  and  vertebral  arteries,  with  arteriolar  branches. 

Capillaries  (C)  -  Represents  the  suit  of  brain  capBIaries,  choroid  plexus  and  arteriolar/venous 
capillaries. 

Venous  Cranium  (Vc  )  -  This  is  a  lake'  of  blood  confined  by  thin  walls  immersed  in  the  brain 
mass.  It  contains  a  controlled  mechanism  to  influence  quick  and  slow  fluid  movement  in  the  brain 
during  stress.  It  comprises  deep  and  superficial  systems,  normally  freely  anastomising. 

Venous  Sinuses  (S)  -  These  are  encased  In  semi-rigid  walls  to  prevent  collapse  in  all  but 
extreme  conditions  of  compression. 

The  Cerebro-Spfnal  Fluid  Is  contained  in  two  compartments: 

Ventricles  (F^ )  -  The  four  ventricles  are  treated  as  a  tree  compartment  with  its  own  inflow 
and  outflow.  Although,  normally,  the  resistance  to  the  latter  is  low  (as  the  aqueduct  of  SiMus 
depends  on  a  puteathre  drive  to  maintain  its  patency),  in  disease,  high  resistance  and  even  occlu¬ 
sion  may  occur  and  Isolate  part  of  the  ventricular  system.  In  such  cases,  we  introduce  in  the 
model  the  extra  ventricular  and  the  spinal  fluid  as  an  additional  compartment. 

Extra  Ventricular  CSF  *  Spinal  Fluid  (F^y )  -  The  cisterns  sulci  and  spinal  fluid  compart¬ 
ment  maintains  a  free  communication  with  the  CSF  extra  ventricular  and,  to  a  lesser  extent,  with 


the  four  main  ventricles. 
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Compartments  representing  the  extra-cranial  venous  (system!-.)  circulation  are: 

Extra-th«cal  venous  plexus  (VgT)  -  A  by-pass  ot  the  systemic  flow  included  in  the  bony  spi¬ 
nal  portion  ot  the  cranial-spinal  compartment.  Normally,  this  compartment  is  much  less  involved 
in  the  circulation  to  the  cerebral  system.  However,  in  special  cases  It  communicates  directly  with 
the  ventricles  through  the  spinal  and  extra  ventricular  fluid.  Therefore,  we  include  It  to  accommo¬ 
date  abnormal  situations. 

Brain  Tissue  (B)  -  This  compartment  is  placed  diagrammatically  as  a  central  one,  emphasiz¬ 
ing  its  vulnerability,  especially  due  to  its  location  between  the  venous  compartments.  All  pres¬ 
sures  acting  on  B  must  be  in  such  balance  as  to  preserve  the  tissue  physically,  yet  allow  and  pro¬ 
mote  optimal  macro  and  micro  circulation  of  fluid  transporting  metabolites. 

The  body  section  of  the  model  is  comprised  of  the  following  compartments: 

Respiratory  System  (R)  -  This  compartment  is  activated  by  the  inhaie/exhale  flux.  Its 
volume  deforms  through  interaction  with  the  heart  system,  abdominal,  superior  vena  cava  and 
body's  arteries. 

Abdominal  (Bp)-  This  compartment  is  subject  to  volume  changes  initiated  by  the  environ¬ 
ment  pressure,  inferior  vena  cava  and  the  respiratory  system. 

The  heart  system  circulation  receives  an  input  perfusion  from  the  superior  vena  cava  and  a 
feedback  input  from  the  Inferior  vena  cava.  The  input  flux  is  injected  to  the  Right  Ventricle  (Rv) 
which  mutually  interacts,  by  volume  deformation,  with  the  respiratory  system,  and  by  flow  with  the 
Pulmonary  (P).  The  pulmonary  communicates  by  mutual  volume  deformation  with  the  respiratory 
compartment  and  (by  flow)  with  the  Left  Ventricle  (lv).  The  left  ventricle  exchanges  volume 
deformations  with  the  respiratory  systems  and  discharges  flow  to  the  body  arteries. 

Body  Arteries  (Ag )  -  This  compartment  receives  Its  inflow  from  the  left  heart  ventricle  and 
introduces  the  major  inflow  to  the  cerebral  arteries.  It  interacts,  via  volume  deformation,  with  the 
respiratory  system,  and  discharges  flow  into  the  superior  and  inferior  vena  cava  compartments. 

Superior  Vena  Cava  (V§p )  -  This  compartment  receives  the  outflow  from  the  cerebral  sys¬ 
tem.  It  also  receives  inflow  from  the  inferior  vena  cava  and  from  body  arteries  and  discharges  it  to 


f 
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the  right  ventricle.  Volume  deformation  interaction  exists  with  the  respiratory  system. 

Interior  Vena  Cava  (V^  )  -  This  compartment  enables  flow  between  body  arteries  and  supe¬ 
rior  vena  cava.  It  interacts  by  volume  deformations  with  the  abdominal,  thus  implicitly  obtaining 
information  from  the  environment  which  is  then  transmitted  as  back  flow  to  the  right  ventricle. 

Next  we  describe  the  resistances  Rtj  (sRji )  and  compliances  (sC^ ),  ascribed  to  the 
various  compartments.  Here,  subscript  ij  in  (  denotes  the  mutual  boundary  of  compartments 
i  and./. 

R \cAa  Carotid  Body  resistance  that  regulates  and  controls  the  inflow  from  the  body  to  the 

cerebral  arteries.  The  control  is  governed  by  the  flux,  QA ,  to  the  cerebral  system. 
CAcb  Compliance  between  the  brain  tissue  and  cranial  arteries  compartments. 

CacFev  Compliance  that  attenuates  the  arterial  pulse  transmitted  by  large  vessels  traversing 
the  basal,  cisternal,  sulcal  spaces. 

R^cC  The  sum  of  the  cranial  arteries,  capillary  and  choroid  plexus  resistances.  The  resis¬ 
tance  associated  with  the  capillary  is  auto-regulatory  control.  This  control  between 
the  arteriolar  and  capillary  vessels,  attenuates  the  systole  artery  pulse. 

Rcvc  Resistance  of  the  Arteriolar/Venous  capillary,  accounting  for  the  pressure  drop 
observed  between  them. 

Rcb  Resistance  of  the  Blood-Brain  barrier  (between  the  capillary  and  the  brain  tissue). 

RcFp/  Endothelial  resistance  of  the  Blood-CSF  barrier.  It  describes  the  choroid  plexus  and 

ependynial  secretion. 

CcFp,  Compliance  manifesting  the  arterial  pulse  transmitted  to  the  ventricle  CSF  by  the 
choroid  plexus  and  the  extra-cellular  fluid. 

JRVc B  Resistance  representing  the  Blood-Brain  barrier  (involved  in  cerebral  oedema). 

Rvcs  A  resistance  comprised  of  two  gates.  First,  the  outlet  from  the  deep  venous  circula¬ 

tion  into  the  straight  sinus  via  the  great  cerebral  vein  of  Galen.  According  to  Le  Gros 


i 


Clark,  an  auto-regulatory  control  may  be  associated  with  this  outflow. 

The  second  component  is  the  superior  cortical  veinous  outflow  into  the  sagittal  sinus. 
An  observable  pressure  drop  exists  across  this  resistance  under  normal  conditions. 

R FnB  Ependumal  resistance  of  the  CSF-Brain  barrier. 

CBF„  Multifrequency  pulsed  compliance,  transmitted  rhythmically  by  the  brain  to  the  ventri¬ 
cle  for  axial  drive  of  CSF. 

RpryFir  Resistance  that  exists  only  in  the  case  of  strong  stenosis,  between  the  ventricles  and 
the  extra  ventricular  compartments. 

RFcvS  Resistance  manifesting  villous  tufts  secretion  into  the  venous  sinus. 

A  low  Compliance  across  the  semi-rigid  sinus  walls. 

RFiv  y  Resistance  that  manifests  the  slow  secretion  area  around  the  spinal  root  sleeves. 

Cjp  v  Compliance  between  the  Intradural  spinal  fluid  and  the  systemic  venous  flux,  across 
the  dura  and  the  extra-thecal  venous  plexus.  This  mediates  postural,  respiratory, 
abdominal  and  other  body  fluctuations,  while  setting  up  the  intra  cranial  pressure 
level. 

RVf-rVsc  Resistance  to  secondary  flow  communication  between  the  exira-thecal  venous 
plexus  and  the  superior  vena  cava. 

Rl/scAg  Resistance  to  inflow  from  body  arteries. 

Cy^jt  Compliance  created  by  the  respiratory  venous  flux.  This  compliance  is  controlled  by 
the  involuntary /voluntary  variations  in  the  respiratory  rhythms. 

RvscRv  Resistance  between  the  superior  vena  cava  and  the  heart  system  via  the  right  ventri¬ 
cle 

RVsc  vic  Back  **°w  resistance  into  the  superior  vena  cava  from  the  inferior  one. 

CrRv’  CrP-  CrLy 

These  are  the  compliances  between  the  respiratory  and  the  heart  system,  comprised 
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of  the  right  ventricle,  pulmonary  circulation  and  the  left  ventricle,  respectively. 

CRAa  Compliance  of  the  respiratory  and  body  arteries. 

Crbd  Compliance  between  the  respiratory  and  the  abdominal.  It  transmits  the  environmen¬ 
tal  effects  to  the  abdominal. 

RrvP  Inflow  resistance  from  the  right  ventricle  to  the  pulmonary  section. 

RpLv  Circulatory  resistance  from  the  pulmonary  to  the  the  right  ventricle. 

Rlvab  Resistance  to  flow  discharged  from  the  left  ventricle  (leaving  the  heart  system)  into 

the  body  arteries. 

RAtVK  Back  flow  resistance,  from  the  body  arteries  into  the  inferior  vena  cava,  to  become  a 
flow  feedback  to  the  heart  system. 

CBdvk  Compliance  between  the  abdominal  and  the  inferior  vena  cava,  transmitting  the 
environmental  effect  to  the  abdominal. 

CgoA  Compliance  due  to  direct  communication  between  the  surrounding  environment  and 
the  abdominal. 

Next,  we  write  the  perfusion  equations  for  the  entire  cerebral-body  oompaitmental  setup. 

3.  FLOW  EQUATIONS 

Let  us  consider  a  single  incompressible  fluid  phase  that  approximately  represents  all  the 
relevant  fluids  in  the  brain-body  system. 

In  writing  the  fluid  flow,  or  balance, equations  we  consider  mass  and  momentum  balances 
for  each  compartment  of  the  cerebral-body  compartment  system. 

The  model  described  here  aocounts  for  an  infinite  environmental  (e.g.,  the  surrounding 
atmosphere,  or  an  experimental  setup  testing  the  influence  of  being  subjected  to  several  g's). 

Essentially,  each  balance  equation  states  that  the  temporal  rate  of  increase  of  either  the  fluid 
mass,  or  its  momentum  in  a  compartment,  is  equal  to  the  amount  of  net  influx  of  that  quantity 
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through  the  compartment's  boundaries  plus  the  external  sources  within  the  compartment.  For  a 
constant  density  fluid,  the  mass  balance  reduces  to  a  volume  balance. 

The  fluid’s  volume  balance  equation  In  compartment  n  takes  the  form 

dV„ 

“j7-  +  L9ni  =Q»  (1) 

* 

where  V„  (t )  denotes  the  volume  d  fluid  In  compartment  n  at  time  t  (equal  to  the 
compartment's  volume),  Qn  Is  the  source  associated  wth  compartment  n ,  and  qni  denotes  the 
flux  flowing  out  of  compartment  n ,  through  the  boundary  ot  the  compartment,  to  its  adjacent  i 
compartment.  This  non  rigidity  of  the  compartment's  boundaries  (-walls)  is  expressed  by  a  com¬ 
pliance  factor,  Cnj .  defined  by 

„  dVn 

nj=d^~)  (2) 

where  pnj  spn  -pj  denotes  the  pressure  difference  between  compartments  n  and  j,  on  both 
sides  of  their  common  wall. 

For  low  Reynolds  number  flow,  the  momentum  balance  of  a  fluid  moving  through  a  capillary 
tube  can  be  shown  to  reduce  to  an  equation  that  expresses  linear  proportionality  between  flux  and 
driving  force.  The  latter  is  composed  of  a  pressure  gradient  and  a  gravity  term.  Here  we  assume 
that  a  similar  expression  governs  the  flux  between  adjacent  compartments.  Hence 

Qni  =  Zm  (Pni  +  Y^m  )  =  (*A  )  I  ni  (3) 

where  (zh )  |  m-  (■  z  |  ^  h  |  m- ),  is  the  conductance  associated  with  the  flow  wall  between 
compartments  n  and  i .  (sHn  -Hi )  is  the  altitude  difference  between  compartments  n  and 

Pni 

i,  Y  is  the  fluid’s  specific  weight,  Aw-  ( =  — —  +  )  denotes  the  piezometric  head  difference 

between  compartments  n  and  i  and  zw-  ( =  fZ^ )  Is  called  the  hydraulic  conductivity  factor 
between  compartments  n  and;. 

Upon  substituting  (2)  and  (3)  into  (1),  we  obtain 
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E  (*A)  I* +X(Cp  )!.,,=«„ 

•  j 


{*) 


"h#f®  (p*~pj)  =P»~Pj)  deno'®s  tf»  temporal  rate  of  increase  In  the  difference  m 

pressure  between  compartments  n  and  Its  adjacent  one.  J. 

Under  certain  abnormal  conditions,  a  number  of  compartments  become  active.  These  are 
the  extra  ventricular  and  spinal  fluid  (Fgy)  and  extra-thecal  venous  plexus  (VCT).  m  the 
mathematical  model,  this  fact  is  expressed  by  IntroArcing  a  parameter  which  can  be  set 
either  to  =  0  for  normal  conditions,  or  to  X**,  =  1  tor  abnormal  ones. 

In  view  of  (4),  we  write  the  loSowing  compart  mental  fluid  balance  equations: 

Ac:  0  =  -(zh)  1^^  +  (*A)U,c  +  (Cp)litcB  +(^Cp)Iac/w 


C :  0  =  -(2A)|XcC  +(zh)\cB  +  (zh)\ CVc  +  (zA )  | CF/y  +  (Cp ) | CFp/ 
vc:  Q  =  ~(zh)\cvc  -(*MIbvc  +  (zA)!vcs  +(Cp)lvcB 
B:  0  =  -(zh )  |  CB  -(zh )  |  FlvB  +(zh )  |  avc+(Cp )  |  BVc+(Cp )  |  bf„+(Cp  )  |  ^ 
S  '  0  =  (ar/i ) | SVac  -  (z/t ) | Vcs  -(**)! FtvS  +  (Cp )  | 

F/v :  0  =  -(zA )  |  cF„+{zh )  I  F„b  HXzh }  |  FnrFtv+(Cp )  |  FiyB  +(Cp )  |  FnrC 

Fev :  0  =  (zA)|y„,s  +  (X*ft)|pwvJC -Q*h)\FryFtY  +  (Cp)\FtvS  + 

+  (XCp )  |  FwvAc  +  (KCp )  |  FfV  VtT 

Vet ’  0  =  (Xzh )  |  v„vK  +  (XCp ) )  v^f^ 

vsc :  0  =  ~(*h )  I  svK  ~  {zh )  |  y^y*.  -  («* )  I  a,  vK  +  (*A )  I  vKxv  - 
+  (XzA  )  I  V*rV*c  -  M* )  I  V«,  +  (Cp  )  I  V„R 

Vic'.  0  =  -(zA)|^V(C+  (zAJIy^  +  (Cp)|V|oBl) 


(5.1) 

(5.2) 

(5.3) 
(54) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 

(5.9) 

(5.10) 
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Ab:  0  =  (*A )  1^^ - (tA)IlvA*  +  (zh)\Atylc  +  (*A)Ia* vK  +  (Cp)\a,r  (S'") 


Ry-  0  =  -(zh )  |  vxSv  +  )lgrp  +  (Cp )  lRt,R  (5.12) 

P-  0  =  ~(z^)\rvp  +  (2^ ) I PLy  +  (Cp)\pr  (S-13) 

Lv:  0  =  -{2/i)|piv  +  (zh)|LvAi  +  (Cp)|iva  (5.14) 

R  :  Qr  =(Cp)\RVac+(Cp)\RjlY+{Cp)\Rp+(Cp)\Riv  +  (Cp)\RBD+(Cp)\RAt  (5.15) 
Bd:  0  =  (Cp ) | Bap  +  (Cp ) | Bd vK+(Cp ) | RdA  .  (5.16) 

Since  the  environment  is  of  infinite  extent,  we  can  write  the  relations 
Qr  ■  CBda  Pbd  (6.1) 

PA  =  0  .  (6.2) 


We  assume  a  Monro-Kellie  postulate  according  to  which  the  (almost)  rigidity  of  the  skull  dic¬ 
tates  that 


(7) 


where  T  is  the  time  period. 

A  stenosis  in  the  passage  between  the  Fjy  and  Fgy  compartments,  initiates  a  build-up  of 
the  compliance  between  the  Fgy  and  Aq  compartments.  This,  In  return,  indicates  that 

XFnF„  =  •  («) 


Under  normal  conditions,  with  free  communication  between  Fjy  and  Fpy  compartments 
b0,f' the  f^aln  ventricle  and  the  extra  ventricular  wffl  merge  into  one  compart¬ 
ment  (theF  compartment),  resulting  in 

hp„  ■  hFgr  m  hF  with  X  =  0  .  (9) 

Not#  that  in  (9)  we  have  assumed  Hpn  =HFgYmHF,  because  only  minor  differences  exist 
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between  their  values. 

Let  us  now  combine  ail  compart  mental  balance  equations  Into  a  global  matrix  form, 
normal  behavior  (l.e.  X*,*  =  0),  and  In  view  ot  (5)  and  (9),  we  obtain 


sf  +|4-fi 


where 


dt  dt  - 


«*->• 


&  =  I*Ac  J*C  .Avc  .Af  .Avjc  <^Vk<^Am  i^Bv  '^iv  >^p  ]? 


0 

0 

0 

0 

0 

0 

cvkr{Cb,-Cb,k) 


Qr 


C\r  _  t» 

-jr-{CB.R-CR)QR 
[Ca*b,(Cb.r~Cr)  +  Ca,r(Cbb~Cbbr)] 


Cr,r 


{Crb-CbbR  )Qr 


(Cs,-Cs.r)Qr 
—=T  (Cr,-Cb,r  Mr 


For  the 


(10) 


(il.i) 


(11-2) 


i  4.  I  ill 

w  P»  V  tfl  K 

N  N  N  N  N  N 


—  N 

N  N 


<o  f** 
N  N 


l  l 


+ 

*  1 

$  i 

i  \ 


h. 

■  V  ^ 

<#  <y  I  g 

i  *,  X  X 

»*  N  o  V 

o  o  o  o 


4,  1.  *,  T  t  i,  I  5  5 


}  §  e  b  $  i  i  1 
ill  l  l  i  t 


i  "> 

f  l 


N  N  N  N 


*■*  •' 

N  N 


lu 


o  o  o  o 


Csjr-Cce  Ci3=-  CDF  Cs.s=Icdf+Cfs+Ccf]  C5.6=-Cre 

C  g  g=Cfs 


-14- 


AiSO, 

Cbd  *  CBbR  +  CbbVk  +  CBdAm  (1 1  -5) 

Cr  u  Crvk  +  Cjtpy  +  Cup  +  +  Crb,,  +  Cjui  (ii.6) 

C  =  CBo  Cr  -  C$oR  .  (i  i  .7) 

Equation  (10)  Is  to  be  solved  for  the  unknown  h  vector,  subject  to  an  initial  condition  of  the 

form 

&(,_0)  =  4o  •  (12) 


Upon  solving  (10),  together  with  (12),  we  obtain  the  solution  for  the  temporal  rates  of  change 
in  piezometric  head,  hBo  and  hR ,  in  the  forms 

*bd  =  1(Cbdr~Cr  )Qr  +(Cr  CBoAt  +CBdR  )hAt  +CR  CBd  v,JlvK+ 

(13) 

+CBDRCRVtchvx:+CBIJR  CRRvhRy+CBDR  CRP  hP  ^BoCRLvhLv  1 

A*  =  7=  1(Cb0-Cbdr  )Qr  +(cbd  Cram+CbdrCbdAb  )^a,+Cbdrcbd  vJivK+ 

(14) 

CB„  CRVtc  ^V^-^Bd^RRv^Rv^Bo  CrP  ^P  +  pB0pRLv^Iw  1  * 

The  compart  mental  balance  equation,  written  in  the  global  compact  form  (10).  involves  con¬ 
ductivities  and  compliances  expressed  by  the  matrices  z  and  C,  respectively.  In  order  to  solve 
for  the  piezometric  head  vector,  4 ,  we  need  to  know  the  z  and  C  values.  We  shad,  next,  discuss 
an  inverse  method  for  estimating  these  parameters. 

4.  PARAMETER  ESTIMATION 

To  predict  the  pressure  and  flux  response  of  the  model  to  external  changes,  the  values  of 
the  parameters  C  and  z  must  be  known.  In  order  to  estimate  them,  we  need  measured  values  of 
pressure  in  the  various  compartments  at  a  sufficient  number  of  points  in  time.  Typical  phasic 


t 

\ 
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pressures  are  depicted  in  Figs.  2-11. 


In  the  present  work,  we  assume  that  all  conductivities  and  compliances  have  the  form  of  a 
time  step  function  i.e.,  they  change  from  one  constant  value  to  another  constant  value,  due  to 
abnormal  situations  such  as  disease.  Because  C  and  z  are  constant  during  long  time  intervals, 
and  k(t)ls  assumed  cyclic,  by  taking  a  temporal  average  of  (10),  i.e.,  integration  over  a  period  ot 
time  divided  by  the  period,  we  obtain 

zh*  =Q*  (15) 

where  h*  and  Q*  denote  mean  values  of  the  piezometric  head  and  of  the  source  term,  respec¬ 
tively.  Since  h*  and  Q*  are  known  averaged  values,  in  order  to  solve  for  z  values,  we  rewrite 
(IS)  in  the  form 

h*zv  =  Q*  (16) 

where  f )  denotes  the  head  difference  between  two  communicating  compartments. 

In  the  cerebral  portion  of  the  model,  we  face  ten  conductivities  against  six  balance  equa¬ 
tions.  However,  by  virtue  of  the  Monro-Kellie  doctrine  which  assumes  absolute  rigidity  of  the  cra¬ 
nial  vault,  actually  only  five  equations  of  the  six  are  independent.  Thus  for  this  redundancy,  five 
additional  conditions  are  needed.  One  of  them  is  the  mean  influx  to  the  CSF  (F)  compartment 
Qp  =  0.3ml/m.n.  (Cutalar,  1968).  Two  more  conditions  are  stipulated  by  physiological  data. 
These  conditions  determine  the  scalar  coefficients  (Sorek  et  al.  1 987)  a  and  P  defined  as 

(a  >  0)  (17) 

(0  £  P  <  ~)  (18) 

where  a  («  10000)  indicates  the  ratio  of  the  vein-venous  sinous  to  the  cerebrospinal  fluid-brain 
barrier,  and  P(  =  1 0-3)  is  the  ratio  of  the  blood-brain  barrier  to  the  cerebrospinal  fluid-brain  bar- 


a  = 


zvcs 

ZFB 


and 


P  = 


ZCB 

ZFB 


-16- 


-17- 


(1-5  )Qi-Q' 

ZAmVK  = - -  •  (20-5) 

AA,V,C 


Hence,  in  view  of  (17)  to  (20),  equation  (16)  is  constructed  of  the  following  matrix  and  vec¬ 
tors 


hBr  “P^CB 

0 

0 

hBVc 

P hcB 

hcvc 

0 

0 

hFB 

0 

hpS 

0 

«*VCS 

~hCVc 

0 

-^BV 

zv  =  I Zfb  •  zCVc •  ZFS  •  ZBvf 


(21.1) 


(21.2) 


e‘=[o,Q*-«;.Q;.o]T 


(21.3) 


Upon  solving  the  inverse  problem  for  the  conductivities,  we  move  to  the  estimation  of  the 
compliances. 

Given  information  of  the  simultaneous  pulse  wave  recordings,  p  ),  and  p  )t  in  the  different 
compartments  and  at  various  times  tk  Jk  =  1 equation  (10)  (or  (5)  and  (6))  now  yields  a 
set  of  k  relations  for  the  compliances,  for  times  (  * 

where  f  )=-—  f). 


(22) 
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£*= 


(zhk ) l^c  -  (2**) \cB-{zhk)\cvc-{zhk)\CF 
(zhk)\cvc  +  (zhk)]BVc  - zhk  |Vcj 
(zA*)|ves  +  (2^*)Ij?s  ~(zA*)lsVj 
(z**)lvcs-<z**)lra-(zA*)  lsvJC+(zA*)|ylcC-(zA*)  \cB-bhk)\cve 
(zhk)\SvK  +  izhk) \VlcVtc  +  (zhk) \AtV$c  -  (zhk) \VscRv 
(zhk)  UjVjc  -  (zhk)  I  VkVsc 
{zhk)  I  VtcRy  -  (2^*  )  I  RyP 
(zhk)\RvP-(zhk)\PLv 
(zhk)\PLv-(zhk)\LvAt 

~{zhk)\A,Vlc-{zhk)  \AaVac  ~(zhk)\AaAc  +  (zhk) \LvAb 

-Qr  +  (2^*  )  I AtVK  -  ( zhk  )  I  VjcVjc 


(23.3) 


Note  that  the  z  values  in  (23.3)  are  already  known. 

v 

By  the  Gauss-Markov  theorem,  the  Q  values  are  derived  as  an  assortment  ot  the  set  of 
informations  through  all  K  time  observations  (Sorek  et  al.  1987b),  namely 

o'-  (nT  it)-1  nT  B  (24) 


where, 


it  = 


/r»i3 


(25.1) 


B  =  lR\R2,..&KhT'K 


(25.2) 
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(26.6) 
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K  . 


£  hVxfio  Kz/i*)Ia,Vk  ~  {Zhk  )  |  V|cVjc) 


'VkBd 


k-1 


£(*$„*„  )2 


(26.7) 


K  . 


£  A*v*  1(zA  * )  I  vkrv  ~  (zh  * )  I  RyP ) 


-*v*  = 


*=1 


£  (Arvr)2 
*=1 


(26.8) 


Cp  = 


£  Apfi [(zA *)lsvp  -  (zA*)| PL  \ 
*=  1  _ _ 

£  (Apfi)2 

*=1 


(26.9) 


if 


£  At«[(2A*)|Ptv-{zA*)|Z/vAii] 

-t  *=1 

’'LyR  ^  '  - 

£  (At*)2 

*=1 


(26.10) 


*  .. 


^AbR 


a£  A^r [(zA * ) | x,vxa-(zft * ) UiAc-(zA* ) |Aj,Vjc-(zA*) U*v,J 


£  (At*)2 


(26.11) 


*=1 


K  .  , 


£  Arbd  +  (2A *  )  I  Aa  vJC  -(zA*)|  v^Vjcl 

nr 

“*Bj)R  — — 

£  (Ars0  )2 


(26.12) 


Thus,  we  have  completed  writing  the  numerical  algorithms  for  estimating  the  z  and  C 
values.  The  mean  and  temporal  values  of  the  fluxes,  pressure  and  pressure  rales  are  derived 
from  phasic  changes  as  depicted  in  Figs.  2-11. 

Note  that  for  the  body  system,  the  compartmental  pressure  and  temporal  rate  of  pressure 
change  ride  on  the  appropriate  respiratory  wave.  However,  for  the  cerebral  system,  this  effect  is 
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attenualed. 

We  assume  a  respiratory  pressure  wave  of  the  form 

Px  =  10  sin  a>xt  (27) 

2ft 

where  <ax  —  is  the  inspiratioiVrespiration  frequency,  and  7V(  =  3.4sec.)  is  taken  as  the 

•Tr 

respiratory  cycle  time. 

We  use  the  following  altitude  values: 


HBd  -110  cm. 

HyK  -  140  cm. 

HVc  -  163  cm. 

Hft  m  135  cm. 

HVk  -  126  cm. 

Hb  ■  163  cm. 

HRv  -  130  cm. 

Hj » 116  cm. 

Hp  *  163  cm. 

HP  m  130  cm. 

Hac  ■  160  cm. 

H§  *  163  cm. 

HLv  “ 130 

He  -  163  cm. 

corresponding  to  a  standing  position.  Also,  we  choose  the  following  values  for  pressure  and  rates 
of  pressure  changes: 

pc  =  30 mm.Hg.  ;  pc  =  0 
Ps  -  2mm.Hg.  ;  p,  =  0 
PvK  =  5mm.Hg.  ;  PvK  =  0  . 

This  concludes  the  calibration  scheme  for  the  brain-body  compartmental  model,  with  con¬ 
ductivities  and  compliances  that  remain  constant  during  long  time  intervals. 

5.  PRESSURE  WAVES 

Once  the  linear  model  described  by  (10)  has  been  calbrated,  it  Is  possible  to  obtain  solu¬ 
tions  for  the  pressure  waves  in  all  brain-body  compartments.  In  incremental  form,  the  solution  of 
the  differential  matrix  equation  (10)  for  h  is 

&(!+*)  =  exp(-T)t&(0-|-1fi(<)]  . 


(28) 


F 


fl 
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Here,  denotes  a  time  dtfterenee  between  the  frontier  time  level,  f m+1, 

and  the  backtime.f"* .  The  matrix  t  reads 

t-Ai  Q-'z  (29) 


The  rational  expansion  of  exp(-r ),  leads  to  the  lotto  wing  formula 


exp(-t ) 


7-<1-0)t 
7  +  0t 


(OS0S1) 


in  which  7  Is  the  unit  matrix.  Substituting  (29)  and  (30)  Into  (28),  yields 


(30) 


&(.-*) * U+*  )_1l7-(1-«)x H&<0-j"1Q<r)]  +|_1fi<0  .  (31) 

The  coefficient  0  controls  the  type  of  numerical  solution  scheme  evolving  in  time.  When 
0  =  0,  we  have  an  explicit  scheme  0  =  1,  an  implcit  scheme,  and  0  <  0  <  1  is  the  mixed 
scheme  case. 


Thus,  with  the  choice  of  0,  the  head  waves  in  the  various  compartments  are  calculated 
according  to  (10).  In  view  of  (3)  and  (31),  we  can  also  evaluate  the  compartmental  pressures 

P(,+*)  =  {(7+0x  rV-d-ejt  +  z-'Q{trH)y  (32) 


This  completes  the  modelling  of  the  brain-body  perfusion  pressure,  as  excited  by  environ¬ 
mental  pressure  and  due  to  Inspiration/expiration  flux. 


6.  CONCLUSION 

A  compartmental  model  tor  perfusion  pressure  Interplay  between  the  cerebral,  respiratory 
and  heart  systems  is  presented.  The  perfusion  pressure  response  is  initiated  by  compartments 
altitude,  environmental  pressure  and  inspiratiotVexpiration  flux. 

Flow  between  compartments  is  governed  by  compartmental  step  function  related  conductivi¬ 
ties  and  compiances  that  take  the  forme  of  time  step  functions. 

The  developed  model  can  be  used  to  guide  the  neurosurgeon  in  the  interpretation  of  possi¬ 
ble  consequences  of  management  methods  apptted  to  abnormal  cases.  Following  are  some 


L 


examples  of  management  maneuvers  that  can  be  interpreted  by  numerical  experiments  of  the 
model. 

(a)  Correction  of  a  fault  If  recognizable  and  accessible,  e.g.,  tumours,  blockage  of  flow  chan¬ 
nels.  This  is  the  ideal  case. 

(b)  Offsetting  of  destructive  tensions  that,  usually,  would  be  short  term  amending,  e.g.,  reducing 
excess  pressure  usually  done  by  shunting  the  CSF  volume. 

(c)  Correcting  mechanisms  that  fail  to  confine  or  maintain  the  system  in  optimal  compliance 
against  repeated  stress  in  the  long  time,  e.g.,  introducing  a  gas  bubble  in  a  sac  into  the  large 
hydrocephalic  ventricle  in  order  to  hold  down  the  pressure  peaks  and  create  better 
Pressure-Volume-Time  relationship. 

In  practice,  the  management  of  (b)  is  on  the  CSF  bulk,  where  as  for  (c)  it  also  involves  the 
blood  volume  as  an  indirect  means  of  promoting  optimum  flux. 

The  model  can  provide  information  for  fault-finding,  target  identification  and  monitoring  the 
results  of  rational  management. 
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Figure  2:  During  one  cycle  of  breathing 
( after  Hamit  et  a l.) 
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TYPICAL  PRESSURE  IN  B 


Figure  3:  During  one  cycle  of  breathing 
( after  Hamit  et  al.) 
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TYPICAL  PRESSURE  IN  Vc 


Figure  4:  During  one  cycle  of  breathing 
( after  Hamit  et  al.) 
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Time  (sec.) 

Figure  6:  During  one  cycle  of  breathing 
( after  Shepherd  et  al.) 
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TYPICAL  PRESSURE  IN  Lv 
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Figure  7:  During  one  cycle  of  breathing 
(after  Wiggers) 
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TYPICAL  PRESSURE  IN  Rv 


Figure  8:  During  one  cycle  of  breathing 
(after  Vernon  et  al.) 
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TYPICAL  Inspiration/ Expiration  FLUX 


Time  (sec.) 


Figure  1 0:  During  one  cycle  of  breathing 
(after  Ben  Haim  et  al.) 
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TYPICAL  AORTIC  FLOW 


Figure  11:  During  one  cycle  of  breathing 
( after  Vernon  et  al.) 
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A  COMPARTMENT AL  BRAIN  MODEL  FOR  CHEMICAL  TRANSPORT 
AND  C02  CONTROLLED  BLOOD  FLOW 

Shaul  Sorek1,  Jacob  Bear2,  Moshe  Feinsod3 

ABSTRACT 

A  compartmental  transport  model  Is  developed,  capable  of  predicting  the  evolution  of  CO2. 
HCOJ  and  H*  in  the  cerebrovascular  system.  In  the  model,  the  transport  ot  these  components  is 
simulated  at  a  subset  of  three  compartments:  CSF,  Capillary-Choroid  Plexus  and  Brain  Tissue, 
belonging  to  a  seven  compartmental  assembly  representing  the  entire  brain.  The  remaining  ones 
are:  Artery,  Vein,  Venous  Sinous  and  Jugular  Bub. 

The  model  accounts  for  advection  associated  with  non  steady  perfusion  fluxes  across  semi- 
pervious  boundaries.  Pressures,  associated  with  perfusion,  are  solved  separately  in  the  seven 
compartment  model.  The  three  compartment  transport  model  also  takes  into  account  changes  in 
compartmental  volume,  due  to  displacement  of  its  boundaries,  diffusion  through  boundaries  and 
rate  of  generation  of  substances  by  chemical  reactions.  A  first  order  reaction  rate  is  assumed  in 
the  Brain  Tissue  compartment.  A  parameter  estimation  method  is  then  developed  to 
boundary  ditfusivities  from  time  averaged  observed  values  of  perfusion  pressure,  tension  of 
carbon  dioxide,  pH  values  and  concentration  of  free  hydrogen  and  bicarbonate  tons.  An  equation 
of  state  for  the  arteries  to  capillaries  influx,  as  a  function  ol  COj  tension  in  the  CSF,  is  then 
suggested  as  a  blood  flow  controller  for  C02  excitations.  Upon  solving  all  coupled  mass  balance 
equations,  and  for  a  pre-evakrated  perfusion  pressures  in  the  artery  and  capillary  co  opart  merits, 
one  can  estimate  the  change  in  arteries  to  capillaries  conductance  at  every  time  step. 

Keywords:  Compartmental  system,  non  steady  perfusion  fluxes,  compfiance,  conductance, 
diffusion,  advection,  chemical  reactions,  molecular  carbon  dioxide,  bicarbonate  ion,  hydrogen  ion, 
parameter  estimation. 
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INTRODUCTION 

Experiments  (Greenberg  el  al.  1978)  show  that  the  flow  of  blood  from  the  arteries  to  the 
capillaries  within  the  cerebrovascular  system,  is  controlled  by  changes  in  the  concentration  of 
carbon  dioxide  (CO2)  within  the  brain  tissue.  In  order  to  investigate  changes  in  this  flow,  in 
response  to  changes  in  CO2  concentration,  a  model  that  simulates  chemical  transport  and  blood 
flow  in  the  brain  system,  was  constructed.  Following  earlier  works  by  the  authors  (Kami  et  al. 
1987,  Sorek  et  al.  I988a-1988c),  a  mu  It  icompart  mental  model  was  selected  tor  this  purpose.  In 
such  a  model,  each  part  of  the  brain,  identified  by  a  certain  function  that  II  fulfills,  is  simulated  as  a 
compartment  characterized  by  unique  values  of  time  dependent  variables,  such  as  pressure  and 
concentration  of  considered  chemical  components.  Adjacent  compartments  interact  with  each 
other,  enabling  flow  of  fluids  (here  blood)  and  chemicals  (e.g.  C02  dissolved  in  the  blood)  to  be 
transported  across  intercomparlmental  boundaries.  Mechanisms  of  species  transport  include 
perfusion  with  the  blood  and  molecular  diffusion.  In  addition,  various  chemical  reactions  involving 
C02.  HCO3  and  FT  take  place,  affecting  their  concentration.  A  model  of  this  kind  is  often  referred 
to  as  a  fully  mixed,  lumped  parameter  model. 

Following  Sorek  et  al.  (1988a)  a  model  composed  of  seven  compartments  was  constructed 
to  simulate  the  following  brain  parts:  Artery  (A),  Capillary  and  Choroid  Plexus  (C).  Venous  (V), 
Venous  Sinus  (S),  Jugular  Bulb  (J),  Cerebrovascular  Fluid  (F)  and  Brain  Tissue  (B). 

In  the  model  investigated  here,  it  is  assumed  that  a  single  (averaged)  incompressible  fluid 
flows  through  the  various  compartments  that  simulate  the  cerebral  system.  The  mathematical 
model  is  comprised  of  mass  and  momentum  balance  equations  for  the  fluid  and  mass  balance 
equations  for  each  of  the  chemical  substances  influencing  the  C02-transport. 

THE  LUMPED  PARAMETER  FLOW  MODEL 

Figure  1  shows  three  adjacent  compartments;  ( i),(j )  and  (n),  separated  by  oommon 
boundaries  which  permit  fluid  to  leak  through  them,  from  (i )  into  (n )  and  from  (n )  into  (J)  They 
also  allow  components  dissolved  in  the  fluid  to  be  transported  across  the  boundaries.  The  solute 
is  transported  through  any  boundary  by  two  mechanisms:  advection  with  the  fluid,  commonly 


referred  to  as  transport  by  perfusion,  and  molecular  diffusion.  The  boundary  itself  acts  as  a  semi- 
pervious  membrane  which  can  be  displaced  in  response  to  pressure  changes  in  the  adjacent 
compartments.  This  displacement  produces  changes  in  the  volume  of  each  compartment.  We 
use  the  term  compliance  to  describe  the  relation  between  the  change  in  a  compartment's  volume 
to  the  change  in  the  pressure  difference  across  It. 

For  a  fluid  of  constant  density,  p,  the  mass  balance  for  the  (n  )th  compartment  is  expressed 
by 

•?*>  t/V), 

TsQin  ~  H  Qnj  +  Q(n)=  j,  (1) 

i-1  j- 1 

where  /(n)  and  </(n)  denote  the  number  of  compartments  from  which  flow  enters  the  (n)th 
compartment,  and  leaves  it,  respectively,  q &  and  qnj  denote  the  fluxes  from  the  (i)th 
compartment  into  the  (n)th  one,  and  from  the  (n)th  one  to  the  (/)t h  one,  respectively,  Q(n) 
denotes  the  fluid  sources  in  the  (n  )th  compartment,  and  Vn  denotes  the  volume  of  compartment 
(n ).  The  compartmental  volume  is  affected  by  the  rigidity  of  the  displaced  boundary  in  response 
to  changes  in  the  pressure  difference ( =  pm  -pn)  across  the  mutual  boundary  between 
compartments  (m )  and  (n ). 

Within  the  cerebral  system,  differences  in  elevation,  H ,  are  small.  Hence,  we  may  assume 
that  between  two  compartments,  m  and  n ,  we  have 

Pmn  =(P«  -Pn)*>Pg{Hm  ~  Hn)  (2) 

where  £  denotes  the  gravity  acceleration. 

The  flux  through  a  semi-pervious  boundary  between  two  compartments  is  proportional  to  the 
piezometric  ( s=  hydraulic)  head  difference  between  them.  However,  in  view  of  (2),  we  assume 
proportionality  to  the  pressure  difference,  i.e. 

9in  =  r&in  Pin  (3) 

where  is  called  the  conductance  of  the  boundary  between  the  (i  )th  and  (n  )th  compartments. 
The  rate  of  change  in  V,,  appearing  in  (1)  is  expressed  by 
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dvn  dVn  dP ^ 

dt  =  Z,  dP«„  '  dt  =^,C-  A 


where  C„,  denotes  the  value  o(  the  compliance  ol  the  (n)th  compartment  due  to  the 
displacement  ot  Hs  boundary  with  the  (m)th  one,  and  is  the  number  of  boundaries  ot 
compartment  (n )  that  are  associated  with  a  compliance. 

The  compliance,  C„„ ,  is  defined  by 


(4.2) 


In  view  of  (3)  and  (4),  we  can  now  rewrite  (1 )  in  the  form 

A»i  dPnm 

Q(n)~~  £  %in  Pin  +  £  %nj Pnj  +  2 

i-1  J-1  m-1  “* 


(5) 


In  general,  the  various  Z's  and  Cnm  s  are  assumed  constant.  In  this  paper  (as  in  earlier 
papers  by  the  authors  (Sorek  et  al.  I988a-1988c),  these  coefficients  are  allowed  to  vary  as 
functions  of  time  e  g.,  in  the  form  of  a  step  function  in  time,  say,  in  response  to  illness.  As  will  be 
shown  in  the  section  below,  the  conductance  between  the  arteries  and  the  capillaries,  depends  on 
the  pressure  differences,  PAC . 

Equation  (5)  constitutes  a  set  of  n  equations,  one  for  each  of  the  n  compartments,  in  the  n 
unknown  pressures,  Pn . 


TRANSPORT  OF  CHEMICAL  COMPONENTS 

Our  next  objective  is  to  formulate  the  mass  balance  equation  for  the  chemical  components 
transported  through  the  compartmental  system. 

We  make  the  assumption  that  the  concentrations  of  the  various  chemical  components  and 
the  various  coefficients  that  govern  transport  and  chemical  reactions  are  independent  of  the 
pressure.  This  makes  the  problem  of  pressure  variations  in  the  compartments  independent  of  the 
problem  of  concentration  distribution.  Accordingly,  the  former  problem  can  be  solved  first,  using 
the  resulting  pressures  as  input  to  the  latter  one.  In  fact,  the  authors  have  already  addressed  the 


problem  of  pressure  distribution  in  earlier  papers  (Sorek  el  at.  (I988a.b)). 

A  detailed  literature  survey  on  trie  effect  of  changes  in  arterial  carbon  dioxide  tension  on 
cerebral  blood  flow  is  documented  in  Greenberg  et  al  (1978).  In  that  work,  a  mathematical  model 
is  developed  that  relates  the  cerebral  blood  flow  control  to  changes  of  CO,  tension  in  the  arteries 
(Pa  CO2).  Greenberg's  model  of  the  controlled  perfusion  in  the  brain  system  is  composed  of  three 
compartments,  capillaries  (blood),  cerebral  extracellular  fluid  and  cerebral  mtracefluiar  fluid.  His 
model  is  based  on  the  assumption  of  rigid  compart  mental  walls  and  constant  perfusions  blood 
pressure. 

In  what  follows,  we  present  a  modification  and  an  extension  of  Greenberg's  model.  The 
extra/intro  cellular  fluid  compartments  are  lumped  into  one  brain  tissue  compartment.  We  relax 
the  assumption  of  rigid  walls  by  introducing  compliances  and  by  solving  the  perfusion  pressure  for 
a  nonsteady  system. 

The  general  layout,  composed  of  a  seven-compartment  system,  is  shown  in  Figure  2.  This 
general  system  was  developed  by  Sorek  et  al.  (1988a)  to  evaluate  changes  in  the  cerebral 
perfusion  pressure.  Figure  3  shows  the  transport  as  conceived  by  Davson  (1967)  and  Greenberg 
et  al.  (1978).  Our  aim  is  to  include  in  the  model  the  regulation  of  arterial  to  capillary  conductance 
as  a  function  of  CO,,  pressure  associated  with  perfusion  flux  and  changes  in  compartment  volume 
related  to  boundary  compliance. 

Following  the  discussion  of  Greenberg  et  at.  (1978),  we  focus  our  attention  on  the  hydration 
reaction  that  takes  place  in  the  physiologic  fluid  circulating  through  the  cerebrovascular  system 

C02  +  H20  H*  +  HCO3  .  (6) 

This  chemical  reaction  influences  the  transport  of  dissolved  C02  and  wilt  thus  be  considered 
in  the  model.  The  coefficients  xy  and  x*  are  the  forward  and  reverse  rate  constants, 
respectively. 

In  view  of  (6),  the  three  chemical  species  that  are  involved  in  controlling  the  blood  flow  are: 
molecular  carbon  dioxide  (C02),  bicarbonate  ion  (HCO^)  and  the  hydrogen  ion  (H*). 
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ln  developing  the  transport  equations,  we  write  a  mass  balance  equation  tor  each  ot  these 
substances,  using  the  notation:  2  =  1  for  CO2,  2  =  2  tor  HCOj  and  2  =  3  for  H*. 


For  the  mixed  celt  idealization  considered  here,  with  no  external  sources  or  sinks  ot  the 
components,  the  mass  balance  tor  the  2  substance  in  compartment  n ,  is  expressed  by 


where  c£  is  the  concentration  i.e.,  mass  of  species  2  per  unit  volume  of  compartment  n ,  Rln 
denotes  the  rate  of  production  of  the  2  th  substance  by  chemical  reaction  in  the  n  th  compartment 
c/(n)  is  the  concentration  of  species  2  in  cell  i  entering  ceil  n ,  Djn  is  the  flux  by  diffusion  of  the 
substance  2  from  compartment  k  into  cell  n  through  their  common  boundary  (a  total  of  such 
boundaries).  The  perfusion  fluxes,  qM  and  qnj.  are  related  to  the  pressure  difference  by  (3).  As 
the  pressure  fluctuates,  these  advective  fluxes,  that  carry  solutes,  also  fluctuate  in  magnitude.  In 
most  existing  models,  these  fluxes  are  assumed  constant.  For  example,  the  flux  from  the  arteries 
to  the  capillaries,  which  affects  the  corresponding  solute  transport,  varies  in  the  range  ±5%  (for  an 
average  of  750  ml/min;  (Sorek  et  al.  f988b)).  The  diffusive  flux  is  expressed  by 


(4(n)-cb  (8) 

where  is  a  coefficient  that  is  related  to  the  diffusivity  and  “thickness"  of  the  boundary,  thus 
regarded  as  boundary  diffusivity.  By  expanding  the  first  term  of  (7)  and  in  view  of  (4),  we  obtain 


f  (c'VB)  = 


dP, 

L  cBM  — ! 

m*1 


dt 


l  f7  dc‘n 

C"+Vn  IT  ■ 


(9) 


By  integrating  (4.2),  we  use  the  symbol  Vn ,  for 

_  o  « 

vn=vn+  I  (Vnm+Cnmpnm) 

m«  1 


(10) 


where 

Vnm=Cnm(h-K) 

and  ,Pn  ,Pm  are  reference  values  given  at  the  same  time  (e  g.,  a lt=0). 
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If  we  now  substitute  (3).  (8)-(l  1 )  into  (7)  we  obtain 


Then,  by  multiplying  (5)  by  cl ,  and  subtracting  from  (12),  we  obtain 


Henry's  law  states  that  the  concentration  ot  a  dissolved  gas  (e.g.  C02)  in  a  fluid  is  related  to 
its  solubility,  of ,  by 

c^=a,4  (14) 

where  of  (constant  for  all  cells)  is  in  moles/liter/mmHg  and  itj,  is  partial  pressure  of  the  gas  in  cell 
n.  In  writing  (14),  we  have  overlooked  the  difference  between  the  volume  ot  a  cell  and  the 
volume  of  fluid  in  it.  assuming  that  the  volume  of  (solid)  tissue  in  the  cell  is  negligible. 


Thus,  in  view  of  (14),  the  mass  balance  equation  (13)  lor  the  dissolved  gas,  becomes 


where  jt/„  =  rt/  -  rt^.  Note  that  (15)  does  not  include  any  eftlux  terms. 

In  view  of  the  transport  processes  depicted  in  Figure  2  and  (15),  we  write  the  explicit  forms 
of  the  mass  balance  equations  for  dissolved  C02  in  the  C,  F  and  B  compartments 


Capillary  and  Choroid  Plexus 

C:  Vc  =  ZAc  *acPac  +  ®bc  kbc  +  ftjlc  +-Rc/a1  (16.1) 


Vc  =  %  +  Vcf  +  Ccf  Pcf 


(16.2) 


where  (VCF  +  CCf  Pcf)!%  5  5%  (Sorek  et  al.  1988a). 
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CSF 

F:  Vf  =  (Zqf Pet  +  )  ^cf  +  ‘Pbf  nBF  +  Pf / ^  (17.1) 

Vjp  =  +  V/sp  +  V)rg  +  Vre  +  CCfPfc  *  CfbP fb  +  CfsPfs  (17.2) 

where  (Vpc  +  Vjb  +  Vfs  +  CcfPfc  +  Cra +  CpsPps )/%  5 1  -5% 

Brain  Tissue 

B :  V"B  -  =  ( ZqbPcb  +  ®bc)  k^b  +  (ZfbPfb  +  ®bf)  ftjla  +  PB  / a  (18.1) 

at 

Vfl  =  ^  +  Vbv  +  Vjjjp  +  CbvPbv  +  CfbPBf  (18  2) 

where  (VBy  +  VBp  +  CBvPbv  +  CFBPbf)/^b  =0.1%.  We  note  the  trend  in  the  effect  of  the 
compliances  on  the  CO,  concentration  changes. 

With  (1 5)  as  a  typical  balance  equation  tor  gas  such  as  C02  in  a  compartment,  expressed  in 
terms  of  the  partial  gas  pressure  it,  we  write  seven  balance  equations  for  the  seven 
compartments.  Compartment  A  serves  as  a  compartment  with  known  fluid  pressure  and  partial 
gas  pressure  that  represent  boundary  conditions.  This  leaves  only  six  compartments  and  six  C02 
balance  equations.  By  examining  Figure  2,  we  note  that  the  three  compartments,  C,  F  and  B. 
contain  no  influx,  excluding  the  flux  from  A  into  C.  Hence,  in  view  of  (15),  we  may  write  a 
subsystem  of  three  balance  equations  for  the  C,  F  and  B  celts  in  terms  of  the  variables  jql ,  k}- 
and  ttB.  This  subset  of  equations  is  independent  of  what  happens  in  the  remaining  three 
compartments  V,  S  and  J.  Actually,  in  order  to  decouple  this  subset  of  equations  from  the 
remaining  equations,  we  had  to  introduce  one  more  simplifying  assumption.  This  assumption 
states  that  the  advective  C02  fluxes  between  the  adjacent  compartments  from  A  to  C,  from  C  to  V, 
and  from  B  to  V,  are  much  larger  than  the  diffusion  fluxes  between  the  same  compartments.  In 
each  case,  we  may  express  this  assumption  by  a  Peclet  number,  P,  =QL  ID ,  that  expresses 
the  ratio  between  the  advective  and  diffusive  fluxes,  where  d  is  the  mean  fluid  velocity  crossing 
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the  appropriate  boundary,  D  is  the  coefficient  of  diffusion  and  L  is  a  characteristic  distance 
between  the  projection  of  compartment  centers  on  the  line  perpendicular  to  the  common 
boundary. 

In  the  compartmental  balance  equations,  the  Pec  let  number  is  expressed  in  the  form 


P.= 


ZinPin 


'Dl 


(19) 


Hence,  whenP,c-*<»,  the  advective  flux  is  dominant,  while  P,c-*0  means  the  domination  ot  the 
diffusive  flux.  When  0  <  Plc  <  we  have  to  include  both  the  advective  and  diffusive  fluxes. 

Our  next  objective  is  to  write  the  balance  equation  for  HCO3  for  the  C,  F  and  B 
compartments.  Here  we  have  to  take  into  account  the  ionic  nature  of  this  component.  In  addition 
to  diffusion  across  the  boundary  between  adjacent  cells,  as  driven  by  a  concentration  difference, 
another  diffusive  flux  is  produced  by  the  difference  in  electric  potential  across  the  boundary. 
These  differences  in  electrical  potential  result  from  the  behavior  of  boundaries  as  semi-pervious 
membranes  which  prevent  certain  ions  from  passing  through  them  (Greenberg  et  al.  1978). 
Because  of  the  electrical  charges,  and  because  we  deal  with  the  diffusive  flux  of  charged  particles 
(viz,  the  HCQj  ion)  the  membrane  itself  becomes  charged  and  affects  the  diffusive  flux  through  it. 

Nernsts  equation  describes  the  equilibrium  potential  created  by  the  ion  distributions  across 
the  membrane. 


The  diffusion  of  an  ion  through  a  charged  membrane  as  governed  by  the  electrochemical 
potential  was  formulated  by  Harris  (Harris,  1960). 

Typical  measured  values  of  the  electrical  potential  across  membranes  for  bicarbonate  and 
hydrogen  ions  (Messeter  and  Siesjo,  1971;  Sorensen,  1971)  do  not  agree  with  values  as 
calculated  by  the  Nemst  equation.  This  leads  to  the  conclusion  that  the  Nemst  equation  does  not 
describe  the  equilibrium  in  electro-chemical  potentials  across  a  membrane  boundary  in  the 
cerebral  fluids,  as  produced  by  the  concentrations  of  HCOi  and  t-F  ions  on  both  sides  of  such 
boundary.  Greenberg  (1978)  addresses  this  Issue  and  suggests  an  additional  diffusive  term, 
controlled  by  a  coefficient  of  diffusion  V" 
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In  what  follows,  we  follow  Greenberg's  (1978)  development,  when  applying  the  component 
balance  equation  (13)  to  HCOf  in  the  C,  F  and  B  compartments. 

Capillary  and  Choroid  Plexus 

C:  Vc  =Zac(c£-c£)Pac +  'Dbc(p$-Y&cc£) 


+  ®Ic(c/_yA:  c£)  +  R£ 


(20.1) 


where 


®ic  =  ®bc  +  ®BC 

yJc  =  ^r- '  exp(-  ~  Em ) 

®BC 


(20.2) 


®FC  =  +  ®FC 

Y&  =  •  exp(- 


KTEm) 


(20.3) 


In  these  equations  /  is  Faraday  constant,  2  is  the  valence  of  ion,  R  is  gas  constant,  T  is 
absolute  temperature  and  Em  is  the  maximum  equilibrium  potential  according  to  the  Nemst 
equation 


Em=  61.5  log 


(21) 


where  c '  and  c  "  are  concentrations  on  both  sides  of  the  membrane. 


CSF 

F :  VF  -  ZCF  (c£  -  c£)Pcf  +  ©Ic  (y!cc(?  -  c/)  + 


+  ©if(ci  -  yIfc/)  +  * / 


(22) 


Here  T>fc  .  Y&:  •  ®  If  arK*Yif  are  similar  in  their  expression  in  (20  2),  (20.3). 
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Brain  Tissue 

»  9 

B :  Vb  =  Zqb  fcc-ciWcB  +  (c/  -  c/JBfb  + 

+  ®Ic(yJcC(?  -  ci)  +  -  c/)  + 


(23) 


Note  that  from  (6),  it  follows  that 

Bn,=-Bn2  =  -B„3  n  =  B  ,C  ,F 


(24) 


provided  that  the  i?'s  are  expressed  in  moles.  Otherwise,  an  appropriate  conversion  must  be 
introduced  (say,  from  moles  to  grams).  Note  that  (24)  constitutes  two  equations  for  each 
compartment.  The  concentration  of  water  is  assumed  constant  so  that  it  has  no  rate  of 
production. 


In  the  case  of  the  Hydrogen  ion,  we  realize  that  a  considerable  quantity  of  the  formed  H* 
becomes  bound  to  various  buffers  (to  hemoglobin  in  the  blood).  The  bound  portion  of  the  total 
amount  does  not  contribute  changes  in  pH  which  is  a  measure  of  the  (non-bound)  concentration 
of  H*  in  the  fluid.  We  define  a  buffer  capacity,  p,  of  a  homogeneous  medium  by  the  rise  in  total 
hydrogen  ion  concentration,c3‘  .per  unit  rise  in  pH,  i.e. 

dc* 


0  = 


dpH 


(25) 


The  definition  of  pH  (  =-logc3),  together  with  (25)  yields 

dc3*  _  -p  dc3 
dt  ~  (In10)c3  dt 


(26) 


The  integration  of  (25)  for  compartments  n  and  i  (assuming  p  to  be  constant  within  the 
range  encountered  in  the  compartments)  yields 

c*  - e*  =  p (pHi  - pHn )  -  (B„  -BJ  (27) 

where  fl„  and  B(  are  integration  constants  associated  with  compartments  n  and  i ,  respectively. 


The  discussion  presented  above  on  the  diffusion  of  the  bicarbonate  ion  across  boundaries 
between  adjacent  compartments  is  also  applicable  to  the  hydrogen  ions. 

In  view  of  (12),  (26)  and  (27),  we  write  the  balance  equations  for  the  H*  ion  in  compartments 
C,  F  and  B. 


Capillary  and  Choroid  Plexus 

C :  VC^=-  — ^  cc3  {  Zac(VpHac  -  Bca  Wac  + 

(28  1) 

+  ®bc  (cb~Sbccc)  +  ®re(c/-8rccc)  -Rc  } 


where 


®Jc  -  ®BC 


/Z 


®bc  exp  (£m  -i-)  +  ®bc 

c3 

8tC  '  ®ic 


(28.2) 


®ir  =  ®rc 

®rc  exp  (£m  +  ®re 

6rc'  sit 


Herep/Zy  -  pH j  and  By  =  Bf  -B^. 


CSF 

dc  ® 

F ;  ^  -rff  =  "  ^  CF  (  ZCF(PPHCF-Brc  )PCF  + 

+  ®fc(5^cC(?-c>?)  +  ®Jf  } 


(28.3) 


(29) 


Here  ®ro-  5^c,  ®Jf,  and  8|f  have  expressions  similar  to  those  given  in  (28.2)  and 
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Brain  Tissua 


B  :  Vg  —Jj-  =  -  lf1p  —  C$  {  Zca  (Pp# CB  -  Bgc  )Pcb  +  %FB  (P pHpB  -BbF  )Pfb  + 

+  ®Jg(SJccC  -  Cfl)  +  (Djjp(8jf'c/-CB“-RB  } 


(30) 


The  presence  of  the  enzyme  carbonic  anhydrase  in  the  red  blood  cells  makes  the 
hydration/dehydration  practically  instantaneous.  Consequently  reaction  (6)  is  always  one  of 
equilibrium,  described  by  the  reaction  constant  k. 

On  the  other  hand,  in  compartment  B,  the  absence  of  carbonic  anhydrase  in  the  cerebral 
extracellular  fluid  means  that  reaction  is  not  one  of  equilibrium.  Instead,  it  has  a  finite  constant 
rate.  We  assume  a  first  order  reaction  in  which  the  rate  is  proportional  to  the  concentration. 
Under  conditions  of  non-equilibrium,  the  rate  of  production  of  C02  equals  the  difference  between 
the  reverse  and  forward  reaction  rates.  Therefore,  in  view  of  (6)  and  (14),  we  write  for  the  B 
compartment 

Rb  -  Vb(kr c$cg  -  tty a’tcjj)  (31) 


In  compartments  C  and  F,  where  the  reaction  is  in  equilibrium,  we  write 


(32) 


Concentration  of  water  (HzO)  is  assumed  constant,  hence  it  is  incorporated  in  the  reaction 
constant  k.  By  virtue  of  (14)  and  (32)  we  conclude  that  c3  =c3(a’rt1  ,c2).  Therefore,  we  can 
replace  (32)  for  compartments  C  and  F,  by  the  relation 

dc3  ot’wt1  dc2  .  a1*  dny 

dt  ~  (c2)2  dt  +  c2  dt  (33) 

At  this  stage  we  have  a  model  consisting  of  18  unknown  variables: 

c„  and  Rl  ;  n  =B  ,C  ,F;  l  =  1,2,3 


for  which  we  have  18  coupled  equations: 
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nine  component  balance  equations,  one  tor  each  component  in  each  compartment, 
two  equations  (33)  for  compartments  C  and  F  and  equation  (31 )  tor  compartment  B,  and 
three  times  equations  (24),  one  for  each  compartment. 

In  principle,  a  solution  can  be  obtained. 

We  note  that  this  set  of  coupled  equations  is  non-linear.  The  coupling  is  through  the  rate  of 
production  terms  (E* ). 

In  the  above  set  of  equations  we  encounter  various  coefficients,  or  parameters.  One  group 
consists  of  coefficients  such  as  try  ,  **,/••  •  thal  are  universal  and  assumed  known.  They  are 
independent  of  the  particular  model  employed  here.  The  second  group  consists  of  the 
conductances  xJm,  compliances  and  boundary  diffusivities  'Dlni  The  determination  of  the 
coefficients  Zni  and  was  described  in  an  earlier  paper  (Sorek  et  al.  I988a,c)  and  will  be 
considered  here  as  already  known. 

In  the  following  section,  we  shall  consider  the  determination  of  the  various  boundary 
diffusivities  of  the  model. 


ESTIMATING  DIFFUSION  COEFFICIENTS 

In  view  of  (16)  to  (18),  we  may  write  the  C02  balance  equations,  in  a  matrix  form 

V  +  ©' it*  =fi’  +S'/a’ 

—  (it  =  — 


(34) 


where 

rt1  =  [*<*,  *ilr 


(35.1) 


(35.2) 
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(%ACPAC  +  ®bc  +  ®rc)  ~®Jc  -®BC 

®'  =  -\Z,cfpcf  +  ®Ac)  (%cfpcf+'Dk:+'£>bf)  -Dbf 

~{ZcbPCB  +  Dbc)  ~(Zpbpfb+'dbf)  (Zcbpcb  +Zfbpfb  +®BC+®Bf ). 

(35.3) 

fi'  =  [-ZUc^AC*A.O,0]r  (35.4) 

S1=lR^.BjJ.«B,Jr  (35.5) 

Similarly,  in  view  of  (20)  to  (23)  we  write  the  balance  equations  for  the  HC03~  element 

V^i  +  V2(L2  =  B2  +  R2  (36) 

—  at  — 


where, 


£2  =  [c<?.  c/.  ]T 

(37.1) 

(ZACP *c+®IcyIc+®Ict/c) 

-®lc 

-‘Die 

C2  =  -{ZcfPcf+Z>fc  y$c) 

(ZCf^>cf+®lc+®l/'Ylf) 

-Dbf 

-{ZcbPcb+'DbciIc) 

-{Zfb  Pfb  +D&f  iIf  ) 

(ZcbP CB+ZFgP fb+Dbc+'Dbf) 

(37.2) 

B2  =  [Zacpacca  •  0,0 f 

(37.3) 

R2  =  [Rt,R?,Ri]T 

(37.4) 

Obviously  the  solution  of  the  model  requires  appropriate  initial  condition.  Examination  of  the 
matrix  forms  (35.3)  and  (37.2)  indicates  that  they  are  non-singular,  which  means  that  a  solution  is 
posstole.  This  is  due  to  the  fact  that  the  model  involves  both  diffusion  and  perfusion. 

The  COj  concentration  (like  those  of  HCO3  and  H*)  continuously  undergoes  changes. 
These  changes  are  introduced  both  by  the  person's  behavior  and  by  various  environmental 


-52- 


f actors,  e.g.,  temperature.  On  the  other  hand,  various  control  mechanisms  exist  in  the  brain  that 
continuously  act  to  restore  some  average  or  normal  concentration. 

As  the  C02  mass  (as  well  as  that  of  HCOs-  and  H*)  within  each  compartment  varies,  we  may 
identify  points  in  time,  say,  t ,  and  1 2  (>f  i ),  at  which  this  mass  is  the  same.  This  means  that  over 
that  time  interval  (regardless  of  the  sequence  of  changes  inbetween  those  points),  <2  -  <1,  the 
total  mass  of  CO2  in  the  compartment  has  not  changed.  Hence 

=0  ;  Vc|,,  =  Vc|r, 

r,  dt 

Thus  by  integrating  (34)  and  (36)  over  a  period  of  time  beginning  and  ending  with  the  same 
component  mass,  we  obtain  accordingly 

©’i’-a’+avo1  os) 


+  (39) 

where  (  ' )  denotes  a  value  averaged  over  the  considered  time  period.  Eq.  (39)  reformulated  to 
solve  the  inverse  problem  for  the  diffusion  parameters  2> 1 ,  take  the  form 

j,®,=L1  (40) 

where 

=  [®ic,  v'fc,  <D'bf)t  (411) 


*CB  *CF  0 
0  it pc  ips 
KjJC  0  kBF 


(41.2) 


,  -1  -  Rc  -1  - 

L'  =  1(Zackacpac  +  ~ r)<  (ZcF^cpPcF  +  ~~t)’ 


(Zcb*CBPCB  +  ZfB*FBPFB  +  ^~)]T 

a 


(41.3) 


where  -  ij)  and  Pi}(  -Pj)  are  differences  between  time  averaged  values  of 


4*. 
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COt  tension  and  pressures  associated  with  perfusion  fluxes,  respectively. 

We  note  that  in  (40).  the  terms  expressing  the  perfusion  fluxes  have  been  incorporated  in 
the  £  matrix,  as  they  are  assumed  known  (from  the  solution  of  the  flow  problem  in  which  all  flow 
parameters  are  known,  Sorek  et  at.  1988a).  However,  from  the  point  of  view  of  diffusion,  the  three 
compartments  act  as  a  closed  system,  the  equations  deserving  diffusion  alone  cannot  be  solved 
unless  we  bring  in  additional  information  ,  e  g.  initial  conditions.  Mathematically  these 
considerations  manifest  themselves  by  the  fact  that  the  matrix  in  (41.2)  is  singular.  As  additional 
information  let  us  assume  that 

©bc  =tf'  (42) 

where  A- 1  is  a  known  value. 


By  virtue  of  (40)  to  (42)  we  obtain 

©rc  =  -  —  IZcfP cp*cf  +  (zcb  Pcb  +K 1  +  zfbPfb*fv  + 

*CF 

~  (RhRl)] 

a ' 


(43) 


1  -  -  R& 

?>BF  =  -  Tj—  KZCbP CB+K')*CB  +  zfbPfb  *fb  +  — f  ] 
It FB  a 


(44) 


Similarly  the  inverse  problem  for  HCOg'  is  solved  by 

C2©2=42  (45) 

where 

©2  =  [©ic-  ©«>  ©iff  (46.1) 


0  -{ffcCc-Cp)  -^b-iIf^f) 
-iy&cCc-tS)  o  -(iIfcf-cb) 

m 


(46.2) 
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L2  =  [(Ztcdchc  +£$)•  {ZCF^CFpCF+Rh> 
{Zcb^CbPcB  +  ZfB^FB^FB  +  ^l)] 

where 


(464) 


To  allow  (he  solution  of  (45),  we  assume,  as  above,  that  the  value  for  ©Jc  is  known. 
'Dbc=K2  (47) 


Hence  in  view  of  (45)  to  (47),  we  obtain 

®fc  =  —f  T~2  VcfPcfCcf  +  ZcbPcb^cb  +  ZfbPfb^fb  + 

Cf-wcCc 

+ ^*(y§c^c-3?)  -  +  -Rj )]  (48) 


®1f  =  -f  2  [ZcbPcbCcB  *  ZfbPfbSfB  +  KHy&cCc-Ci)  +  (49) 

The  non-linear  equations  (28)  to  (30)  describe  the  mass  balance  ot  H+(  /  *3)  in  the  C,  F  and 
B  compartments. 

Again,  we  assume  that  the  mass  of  the  hydrogen  ion  undergoes  continuous  changes.  The 
integration  of  (28)  to  (30)  over  a  period  al  the  end  of  which  the  mass  returns  to  its  initial  value,  will 
eliminate  the  time  derivative  and  yield  a  linear  set  of  equations  that  take  the  form 

C3V3  =  L3  (50) 

where 


©3  =  [©lc.©#c.©lFf 


(51.1) 


-{e^ej-Sic  fa?)*!  -(e£ci?-®fc(c£)*l 


0 


g3  *  0  48re  (c£)!  ~  ~{c£c£  -  &lrc3c£\  (51.2) 

^2c(c<?)2-c<?c|]  0  -{Si/-C(?c/  -  c^ca3] 

£t3  =  [2ac(8p#ac  +  Bac)Pacc£  +  Rccc)  • 

.  +  Bcf)Pcfc?  +  Rpc £] ,  (51 .3) 

.  [2cb(Pp^cb  +  Bcb)Pcbc£  +  2/8  (Pp-Hra  +  Bfb)Pfbc£  +  Rbc$]  J7 
To  solve  for  ©3  in  (50),  we  assume,  as  above,  that 

®bg  =  *3  (52) 

where  if3  is  a  known  value. 


By  virtue  of  (50)  to  (52)  we  obtain 


®PC  =  .  - ~  - TT  +  2Cb(PP#CB+BCB  WcB^i  ' 

cSc?-&tc{c$f 

+  Zra  (Ppi/ra  +flra )PracB3  -  K3(cgci-b$c{cSf)  -  (Rfc$+Rbc£)] 


©!jp  = 


lZCB  (VplicB  +B  CB  W CB  c£  + 


cc3cb3-«IbCc3c/ 

+  %FB  (VpHfB+Bf 8  WfBc£  +  Ri(^c(^c)Z~c£c£)  +  #BC^1 


After  estimating  the  diffusion  parameters,  one  still  has  to  consider  additional  parameters. 
These  are  shown  in  Table  1. 

With  the  known  parameters  we  solve  for  the  transport  of  C02  ,  HC03‘  and  H+  in  the 
Capillaries,  CSF  and  Brain  Tissue  compartments.  The  transport  and  flow  problems  are  solved 
simultaneously. 
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BLOOD  FLOW  CONTROL 

In  writing  the  balance  equation  (5),  it  is  usually  assumed  that  the  conductances,  Z,  are 
independent  of  time.  However,  in  reaBty,  the  values  of  Z  vary  continuously  in  response  to  certain 
changes  that  occur  in  the  cerebrovascular  system.  In  what  follows  the  brain  cells  are  continuously 
metabolizing,  consuming  oxygen  and  producing  carbon  dioxide.  The  changes  in  the  C02 
production  are  probably  the  dominant  factor  in  cerebral  blood  flow  control. 

We  use  a  relation  (Ponten  and  Siergo,  1966)  between  C02  tension  in  the  CSF  (it/),  and  in 
the  arteries  (*x). 

ic/  =  8.6  +  0.942  it,}  (55) 


Following  Greenberg  et  al.  (1978),  we  also  use  the  relation  between  arterial  to  capillary 
blood  flow,  and  it,}  (Reivich,  1964). 


Qac  —  ZacPac  =  20.9  + 


92.8 

1+10570  («a)2281 


(56) 


An  increase  in  the  amount  of  C02  in  the  arterial  blood  bathing  the  brain  causes  a  dilation  of 
blood  vessels  and  an  increase  in  the  blood  flow  which  carries  away  C02  and  reestablishes 
normal  conditions. 


By  combining  (55)  and  (56)  we  find  the  relation  between  it/  and  ZACPAC  ■  Thus  for  any 
time  instance,  knowing  the  solution  for  the  perfusion  pressure  fall,  Paq  (i.e.,  by  solving  the  fluid's 
flow  equations)  and  COg  tension  in  the  CSF,  it/,  (i.e.,  by  solving  the  solute  transport  problem),  we 
find  the  appropriate  arterial-capillary  conductance 


Zac  = 


Pac 


20.9  + 


92.8 


1  +  10570  ( 


*—2.281 


0.942 


(57) 


Hence  we  find  the  change  of  arterial-capillary  conductance  due  to  changes  in  perfusion 
pressure  and  C02  tension. 


i 

i 
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SUMMARY  AND  CONCLUSIONS 

A  compart  mental  brain  model  was  developed  to  describe  transport  processes  of  molecular 
carbon  dioxide  (CO2),  bicarbonate  ion  (HCOj)  and  hydrogen  ion  (H*)  between  Brain  Tissue  (B). 
Capillaries  -  Choroid  Plexus  (C)  and  the  CSF  (F)  compartments.  The  model  consists  of  mass 
balance  equations  (i.e.  transport  equations)  accounting  for  non-steady  advection,  diffusion  and 
generation  of  constituents  due  to  chemical  reactions. 

Balance  equations  are  expressed  in  tension  values,  of  CO?,  concentration  for  HCO3,  and  pH 
values  and  free  concentration  for  H\ 

The  fluid  is  assumed  to  be  of  a  constant  density,  hence  allowing  decoupling  between  the 
solutions  of  fluid  and  of  component  balance  equations.  A  separate  solution  is  thus  obtained  for 
non  steady  flow  through  the  cerebrovascular  system  model  by  seven  compartments:  Artery, 
Capillary  plus  Choroid  Plexus,  Brain  Tissue,  CSF,  Vein,  Venous  Sinous  and  Jugular  Bulb.  This 
solution  describes  the  evolution  of  compartmental  perfusion  pressure  subject  to  conductances 
expressing  the  case  of  leakage  through  semi-pervious  boundaries  and  compliances  expressing 
compartmental  volume  rate  of  change  due  to  boundary  displacements  (Sorek  et  al. 
I988a-I988c).  The  non  steady  flow  governs  the  advection  term  in  the  model  describing  the 
transport  of  chemical  components.  This  enables  the  study  of  pathological  situations  such  as 
different  patterns  of  initial  fluxes  that  give  rise  to  transient  periods  of  transport:  and/or  changes  in 
conductances  and  compliances  because  of  occlusions. 

Generation  of  components  in  the  physiological  fluid  within  a  compartment,  is  governed  by 
the  stochiometric  relation  expressing  the  hydration  reaction.  Within  the  capillaries  -  choroid 
plexus  and  the  CSF  compartments,  the  reaction  is  considered  to  be  instantaneous,  hence  in 
equilibrium.  Within  the  brain  tissue  compartment,  the  production  of  CO?  is  assumed  to  equal  the 
difference  between  the  reverse  and  forward  reaction  rates. 

In  order  to  solve  for  the  different  components,  we  need  to  know  various  coefficients.  Some 
of  them  may  be  found  in  different  citations.  However,  most  of  the  boundary  diffusivities  need  still 
to  be  estimated.  Therefore,  we  present  a  parameter  estimation  method  to  assess  time  averaged 
boundary  diffusivities  based  on  time  observations  of  perfusion  pressures,  CO?  tension,  HCOg 
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concentration  and  pH  and  tree  H*  concentration.  For  future  applications,  an  assessment  is  still 
needed  of  the  sensitivity  of  the  transport  solution  to  parameters  such  as  diffusivities, 
conductances  and  compliances  defined  at  the  co mpart mental  boundaries.  For  example,  based  on 
previous  work  of  the  authors  (Sorek  et  al.  I988a-I988c),  we  note  that  the  relation  between 
temporal  volume  change  (expressed  by  the  product  of  boundary  compliance  and  perfusion 
pressure  drop  -  Cy  P,j)  and  average  boundary  volume,  is  much  more  significant  in  the 
capillaries  than  in  the  brain  tissue. 

Diffusivities,  compliances  and  conductances  are  assumed  as  step  functions  in  time  i.e.  they 
may  change  from  one  constant  value  to  another  because  of  reasons  such  as  illness,  aging,  etc. 
We  relax  this  assumption  for  the  conductance,  ZAC ,  controlling  flow  from  the  arteries  to  the 
capillaries.  We  introduce  an  equation  of  state  that  relates  this  conductance  to  perfusion  pressure 
drop  between  arteries  and  capillaries  and  C02  tension  in  the  CSF.  Hence,  ZAC  is  changing 
continuously  in  response  to  excitations  from  the  perfusion  pressure  and  C02  tension.  This  short 
time  of  response  of  ZAC  describes  flow  control  from  arteries  to  capillaries. 
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Table  1:  (Continued) 


Physical  and  chemical  parameters 


Buffering  parameters 
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Figure  2: 

Cerebral  seven  compartments l  system  for  perfusion  pressure  and  transport  processes. 


Schematic  of  the  brain  model  compartments  showing  both  fluid  flow  (solid  arrows)  and 
diffusion  (broken  arrows)  between  compartments.  Transport  Is  (a)  between  extracellular 
fluid  (ECF)  and  capillary  blood,  (b)  between  the  capillary  blood  and  cerebrospinal  fluid 
(CSF),  (c.d)  between  the  CSF  and  ECF,  (e,f)  between  the  cells  (glia  and  neurons)  and  the 
ECF  (modified  from  Davson,  1967 and  after  Greenberg  et  al.  1978). 
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